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SUMMARY  AND  CONCLUSIONS 


This  study  has  produced  a  design  equation  that 
can  be  used  to  estimate  nitrogen  removal  in  waste- 
water  pond  systems.  The  equation  was  developed 
and  validated  with  data  from  full-scale  systems  in 
various  parts  of  the  United  States  and  Canada. 
The  ranges  of  the  parameters  included  in  the  de¬ 
velopment  and  validation  were: 


Wastewater  nitrogen 
Detention  time 
Pond  temperature 
Pond  pH 


8  to  61  mg/L 
5  to  231  days 
l°to  28  °C 
6.4  to  9.5 


As  a  result  the  design  equation  should  be  valid  for 
wastewater  ponds  located  anywhere  in  the  temper¬ 
ate,  subarctic  or  arctic  zones. 

The  equation  can  be  used  in  those  situations 
where  nitrogen  removal  and/or  conversion  is  re¬ 
quired  prior  to  discharge.  It  will  produce  an  esti¬ 
mate  of  nitrogen  removal  that  will  occur  in  a  pond 
system  for  cost  effectiveness  comparison  with  al¬ 
ternative  removal  methods. 

The  equation  should  also  be  used  in  the  design 
of  all  land  treatment  systems  that  incorporate  a 
pond  for  treatment  or  storage.  Nitrogen  is  typi¬ 
cally  the  controlling  parameter  in  system  size  and 
therefore  cost,  and  a  reliable  estimate  of  nitrogen 
removal  in  the  pond  will  allow  more  cost-effective 
design  of  the  land  treatment  component. 

As  a  result  of  this  study  it  is  concluded  that: 

1.  The  nitrate  concentration  in  the  effluent  from 
non-aerated  and  partial  mix  aerated  ponds  will  be 
low.  Essentially  all  of  the  nitrogen  in  the  effluent 
will  be  in  either  the  ammonia  or  organic  form. 

2.  The  nitrogen  removal  in  these  pond  systems  is 
dependent  on  pH,  temperature  and  detention 
time.  Under  warm,  sunny  conditions  with  ade¬ 
quate  detention  and  alkalinity,  nitrogen  removal 


can  exceed  80%.  Although  the  remaining  nitrogen 
will  not  be  in  the  nitrate  form,  it  may  still  be  possi¬ 
ble  to  meet  discharge  standards  because  of  the 
high  overall  removal. 

3.  The  dependence  on  pH,  temperature  and  de¬ 
tention  time  suggests  that  volatilization  (loss  of 
gaseous  ammonia  to  the  atmosphere)  is  the  major 
factor  responsible  for  long-term  permanent  nitro¬ 
gen  removal  from  pond  systems.  However,  there 
is  no  direct,  absolute  proof  to  support  this  hypo¬ 
thesis. 

4.  The  loss  of  gaseous  ammonia  to  the  atmo¬ 
sphere  from  beneath  a  continuous  ice  cover  seems 
unlikely.  However,  nitrogen  removal  was  ob¬ 
served  under  such  conditions  in  the  systems  in  the 
northern  United  States  and  Canada  that  were  in¬ 
cluded  in  this  study.  This  suggests  that  some  other 
mechanism  may  prevail  in  the  winter,  with  volatil¬ 
ization  again  occurring  when  conditions  are  favor¬ 
able. 

5.  The  presence  of  an  ice  cover  on  a  wastewater 
pond  restricts  oxygen  transfer  and  results  in  an  in¬ 
crease  of  dissolved  constituents  in  the  remaining 
liquid  due  to  rejection  during  the  freezing  process. 
An  equation  to  estimate  the  effect  was  developed 
in  this  study.  All  of  these  factors  combined  mean 
that  ice-covered,  non-aerated  ponds  will  produce 
very  poor  effluent  quality  during  the  period  of  ice 
cover. 

6.  Methods  for  estimating  the  pH  and  the  temp¬ 
erature  in  pond  systems  are  presented  in  this  re¬ 
port.  Such  data  are  essential  to  determining  nitro¬ 
gen  removal  and  are  useful  for  other  purposes  as 
well. 

7.  Use  of  the  design  equation  should  ensure 
maximum  cost  effectiveness  for  design  of  dis¬ 
charging  systems  that  require  nitrogen  removal  or 
conversion,  and  for  all  land  treatment  designs 
where  nitrogen  is  the  controlling  parameter. 


NITROGEN  REMOVAL  IN 
WASTEWATER  PONDS 


Sherwood  C.  Reed 


INTRODUCTION 

Wastewater  ponds  are  a  relatively  common  and 
cost-effective  treatment  method.  Labor  and  ener¬ 
gy  requirements  are  low  compared  to  other  me¬ 
chanical  treatment  systems,  and  the  concept  is 
particularly  attractive  for  military  facilities,  where 
sufficient  land  is  usually  available. 

Ponds  can  be  used  as  the  sole  treatment  method 
with  discharge  to  a  nearby  watercourse,  as  a  final 
polishing  step  after  other  forms  of  treatment,  or 
as  preliminary  treatment/storage  prior  to  land  ap¬ 
plication  of  wastewater.  Ponds  are  also  classified 
with  respect  to  the  degree  of  mechanical  mixing  or 
aeration  provided.  Table  1  summarizes  the  charac¬ 
teristics  of  the  pond  systems  in  common  use. 

The  forms  and  transformations  of  wastewater 
nitrogen  are  critical  to  engineering  design  and 
water  quality  decisions.  Nitrogen  in  any  form  may 
eventually  promote  eutrophication  in  receiving 
waters,  and  in  some  cases  its  significant  reduction 
or  total  removal  may  be  required  prior  to  dis¬ 
charge.  In  many  situations,  the  presence  of  am¬ 
monia  or  ammonium  ions  is  critical,  since  their 


oxidation  will  reduce  the  oxygen  content  of  the  re¬ 
ceiving  waters.  In  addition,  ammonia  is  toxic  to  a 
wide  variety  of  fish  in  relatively  low  concentra¬ 
tion,”  and  the  U.S.  Environmental  Protection 
Agency  recommends  a  limit  of  0.02  mg/L  based 
on  its  impact  on  salmonid  fish  fry.  As  a  result, 
many  states  limit  the  discharge  of  ammonia  nitro¬ 
gen  at  least  during  the  spring  and  summer  months, 
when  removal  is  required  prior  to  discharge.  Ni¬ 
trogen  is  also  the  limiting  design  parameter  for 
many  land  treatment  systems;  therefore,  the  con¬ 
centration  of  nitrogen  leaving  the  preliminary 
treatment/storage  pond  is  a  critical  factor  in  de¬ 
termining  the  size  and  therefore  the  cost  of  the 
overall  system. 

It  is  still  a  common  engineering  practice  to  ig¬ 
nore  any  nitrogen  transformations  or  removal  that 
might  occur  in  a  pond  system.  It  has  been  ob¬ 
served  on  numerous  occasions  that  nitrogen  losses 
do  occur  in  pond  systems,  but  there  were  insuffi¬ 
cient  quantitative  data  to  permit  the  development 
of  a  rational  design  approach.  The  recent  publica¬ 
tion"0  of  very  detailed  pond  studies  docu¬ 
ments  nitrogen  removal  of  over  80%  under  f^>or- 


Table  1.  Design  factors  for  treatment  ponds 


Type 

Detention  time 
(days) 

Depth 

(m) 

BOD  loading 
(kg/ha  day) 

Oxidation 

10-40 

0.45-0.91 

67-134 

Facultative 

Winter  average  air  temperature 

above  60  °F 

25-40 

0.91-1.5 

45-90 

32°F-60°F 

40-60 

1.2-1. 8 

22-45 

below  32  °F 

80-180 

1. 5-2.1 

11-22 

Partial  mix  aerated 

7-20 

2. 4-3.0 

34-112 

able  conditions,  and  there  is  now  a  sufficient  body 
of  data  to  develop  and  verify  a  mathematical  de¬ 
scription.  This  report  presents  the  development  of 
a  rational  design  approach.  Use  of  this  approach 
will  permit  the  estimation  during  design  of  the  ni¬ 
trogen  concentration  to  be  expected  in  a  pond  ef¬ 
fluent,  and  thereby  ensure  cost  effectiveness  in  de¬ 
sign,  construction  and  operation  of  any  subse¬ 
quent  treatments  that  may  be  included. 


BACKGROUND 

The  loss  of  ammonia  nitrogen  from  streams  and 
from  alkaline  water  impoundments  was  described 
by  Stratton  in  1968-69." 24  He  suggested  a  first- 
order  reaction  to  describe  the  loss  and  indicated 
that  the  rate  constant  was  proportional  to  pH, 
temperature,  detention  time  and  the  degree  of  air 
movement  and  water  turbulence  at  the  water  sur¬ 
face.  He  estimated  that  the  ammonia  losses  from  a 
treated  wastewater  impoundment  in  Santee,  Cali¬ 
fornia,  approached  6%  per  day.24 

In  a  1968  summary  of  nitrogen  removal 
methods,  Johnson12  suggested  that  the  removal  of 
nitrogen  in  wastewater  ponds  is  dependent  on  the 
growth  and  harvest  of  algae  and  on  soil  adsorp¬ 
tion  of  ammonium  ions  if  seepage  is  allowed.  Al¬ 
gal  harvesting  was  considered  essential  by  the  in¬ 
vestigators  he  cited.  Otherwise  the  algae  would  die 
and  decompose,  and  the  nitrogen  would  be  re¬ 
cycled  to  the  pond  or  to  the  receiving  stream.  Pilot 
studies  by  DiGiano  indicated  that  algal  uptake, 
not  nitrification,  was  the  major  factor  for  the  ni¬ 
trogen  transformations  he  observed.1  Porcella  et 
a!.,1’  in  a  study  of  a  tertiary  effluent  impound¬ 
ment  in  Alpine  County,  California,  noted  that 
48%  of  the  influent  nitrogen  over  a  2-year  period 
could  not  be  accounted  for,  and  they  suggested  ni¬ 
trification-denitrification  as  the  most  plausible 
reason.  They  also  noted  algae  growth  and  sug¬ 
gested  that  trout  mortality  was  due  to  ammonia 
toxicity. 

Ammonia  removals  of  about  70%  were  report¬ 
ed  by  Mara  et  al.14  for  a  five-cell  pond  in  Brazil 
(average  temperature  25  °C,  median  pH  8,  deten¬ 
tion  time  29  days),  but  removal  mechanisms  were 
not  discussed.  Very  significant  nitrogen  removals 
were  measured  in  the  four  EPA  studies  previously 
mentioned'  14  21  but  no  discussion  of  removal 
mechanisms  was  presented.  Ferrara  and  Harle- 
man7  developed  a  mathematical  model  (complete 
mix — dynamic)  that  indicated  that  the  major  path¬ 
way  for  nitrogen  removal  in  facultative  ponds  was 
biological  and  algal  uptake  followed  by  loss  to 


bottom  sediments.  King11  reported  on  nitrogen  re¬ 
moval  in  a  pond  system  in  Michigan,  and  indi¬ 
cated  that  the  harvest  of  aquatic  plants  (algae, 
etc.)  could  only  account  for  about  10%  of  the  ni¬ 
trogen  removed.  He  achieved  up  to  95%  nitrogen 
removal  in  the  120-day  detention  time,  and  sug¬ 
gested  that  loss  of  ammonia,  as  a  gas,  to  th»  at¬ 
mosphere  was  the  major  factor.  He  presented  the 
following  first-order  equation  to  describe  the  re¬ 
sults  achieved: 

N,  =  N0e-°03' 

where  N,  =  nitrogen  (total,  as  N)  at  time  t 
(mg/L) 

N0  =  nitrogen  (total,  as  N)  entering  pond 
(mg/L) 

/  =  detention  time  (days). 

Idelovitch"  described  very  high  ammonia  re¬ 
movals  from  ponds  in  Israel  during  the  period 
1975-79.  The  pH  was  first  raised,  by  lime  treat¬ 
ment,  to  about  1 1,  and  ammonia  removals  of  95% 
were  observed  after  14  days  detention,  with  air 
stripping  identified  as  the  responsible  mechanism. 
In  addition  to  the  alkalinity/pH  requirements, 
Idelovitch  indicated  that  removal  efficiency  is  also 
dependent  on  water  temperature  and  detention 
time. 

It  is  clear  from  the  review  presented  above  that 
nitrogen  losses  from  pond  systems  have  been  re¬ 
cognized  for  a  number  of  years,  but  that  there  has 
been  no  consensus  on  identification  of  the  re¬ 
sponsible  mechanisms.  The  publication  of  exten¬ 
sive  new  data1 10 1621  from  actual  pond  systems  has 
not  yet  resolved  the  issue.  Pano  and  Middle- 
brooks'*  developed  a  complete-mix  mathematical 
model  using  data  from  the  first  cell  in  three  of  the 
EPA  pond  studies' 1621  and  then  applied  the  model 
to  the  results  from  cells  2-4  in  these  same  systems 
for  verification.  Independent  data  from  other  full- 
scale  operating  systems  were  also  found  to  fit  the 
model.  Their  relationship  is  dependent  on  pH, 
temperature  and  hydraulic  loading  rate,  and  air 
stripping  of  ammonia  is  suggested  as  the  major  re¬ 
moval  pathway.  Ferrara  and  Avci6  applied  a 
model  previously  developed  at  MIT7  to  some  of 
the  same  data21  used  by  Middlebrooks  and  Pano, 
but  concluded  that  ammonia  volatilization  (air 
stripping)  accounted  for  a  very  small  fraction  of 
the  total  nitrogen  removed.  They  suggested  bio¬ 
logical  and  algal  uptake  followed  by  loss  as  sedi¬ 
ment  as  the  major  pathways. 

This  continuing  controversy  may  not  be  capable 
of  resolution  with  the  present  data  base.  The  tern- 


perature,  pH  and  detention  time  conditions  that 
might  result  in  high  rates  of  ammonia  volatiliza¬ 
tion  would  also  support  high  rates  of  biological 
activity  and  algal  uptake  of  ammonia.  It  may  not 
be  possible  to  determine  whether  the  Pano  and 
Middlebrooks  model  or  the  Ferrara  and  Avci 
model  is  the  more  valid  description  of  nitrogen  re¬ 
moval.  Additional  concerns  arise  since  both  are 
based  on  a  complete-mix  assumption,  which  is  not 
entirely  valid  for  the  long  detention  time, 
multiple-cell  facultative  ponds  considered.  In  ad¬ 
dition,  the  Pano  and  Middlebrooks  equation  deals 
only  with  the  ammonia  fraction  entering  a  pond 
and  ignores  the  other  nitrogen  forms  that  may  be 
present.  These  other  forms  may  then  cycle 
through  several  stages  and  appear  as  ammonia 
while  the  wastewater  is  in  the  pond.  The  Ferrara 
and  Avci  equations  assume  a  very  significant  ni¬ 
trogen  loss  to  bottom  sediments.  Assumptions  re¬ 
garding  the  permanence  of  this  disposition  seem 
optimistic  and  tend  to  minimize  the  solubility  of 
nitrogen  and  to  ignore  the  resuspension  of  particu¬ 
late  matter  caused  by  density  currents  and/or  sea¬ 
sonal  “overturns”  due  to  thermal  effects  in  a 
pond. 


THEORY 
1.  Nitrogen  forms 

Nitrogen  can  exist  in  wastewaters  in  a  number 
of  different  forms  because  of  the  various  oxi¬ 
dation  states  represented,  and  it  can  readily  change 
from  one  state  to  another,  depending  on  the  oxi¬ 
dation  or  reduction  reactions  induced  by  the 
physical  or  biochemical  conditions  present  at  the 
time.  The  principal  forms  of  concern  in  waste- 
water  are: 

Organic  nitrogen:  Organic  nitrogen  is  present 
in  the  protein  molecules  in  animal  tissue  and  in 
fecal  matter  (typical  assumption  is  about  8-10% 
nitrogen  content).  Microorganisms  decompose  the 
organic  nitrogen  and  produce  ammonia  (NHj- 
NH.  +  ). 

Ammonia:  Ammonia  can  be  present  as  mole¬ 
cular  ammonia  (NHj)  or  as  ammonium  ions 
(NH,  +  ).  The  nitrogen  in  urine  is  initially  urea, 
which  is  hydrolyzed  by  an  enzyme  to  ammonium 
carbonate  ((NH.)iCOj).  Microorganisms  (Nitro- 
somonas  and  Nitrobacter)  oxidize  these  ammonia 
compounds  to  nitrite  (an  unstable  intermediate) 
and  then  nitrate.  Ammonia  can  also  be  taken  up 
by  either  bacteria  or  plants  and  converted  to  cell 
matter,  which  returns  the  nitrogen  to  the  organic 
form. 


Nitrate:  Nitrate  (NOr)  is  the  stable  end  point 
of  the  oxidation  reactions: 

Organic  N  —Ammonia-Nitrite— Nitrate 

Nitrate  can  be  taken  up  by  plants  and  algae  and 
converted  to  organic  nitrogen  in  cell  tissue.  Or,  in 
an  oxygen-deficient  or  anaerobic  environment 
(with  sufficient  carbon  present),  it  can  be  reduced 
to  elemental  nitrogen  (N2)  and  lost  as  gas.  The  ni¬ 
trifying  organisms  do  not  compete  and  perform 
effectively  until  the  simple  carbonaceous  organic 
materials  (BOD)  have  been  substantially  removed 
from  the  wastewater.  It  is  also  believed  that  the  ni¬ 
trifying  organisms  function  more  effectively  if 
they  are  attached  to  surfaces  rather  than  sus¬ 
pended  (e.g.  free-floating)  in  the  liquid  column. 
As  a  result,  a  relatively  low  degree  of  nitrate  con¬ 
version  seems  to  occur  in  typical  wastewater  treat¬ 
ment  systems.  To  achieve  nitrification  the  system 
design  must  provide  for  either  relatively  low  initial 
organic  loadings  or  long  detention  times,  suffi¬ 
cient  oxygen  and  suitable  temperatures,  and,  for 
maximum  conversion  efficiency,  a  surface  for  at¬ 
tached  growth.  As  discussed  in  more  detail  later, 
these  conditions  are  not  simultaneously  present  in 
the  typical  facultative  pond  and  the  apparent  rates 
of  nitrate  conversion  are  low. 

The  nitrogen  concentration  in  typical  municipal 
wastewaters  ranges  from  about  15  to  over  50 
mg/L.  About  60%  of  this  is  in  the  ammonia  form, 
about  40%  in  the  organic  form.25  To  avoid  con¬ 
fusion  it  is  conventional  to  express  these  concen¬ 
trations  as  equivalent  elemental  nitrogen  (as  N)  re¬ 
gardless  of  the  particular  ionic  form.  For  example, 
10  mg/L  of  nitrate  nitrogen  (as  N)  is  equal  to 
about  45  mg/L  of  the  specific  nitrate  compound 
(as  NOr). 

The  presence  of  ammonia  gas  (NHj)  or  ammon¬ 
ium  ions  (NH« +)  in  the  wastewater  is  strongly  de¬ 
pendent  on  the  pH  of  the  water  and  on  tempera¬ 
ture.  Figure  I  illustrates  this  relationship.  At  pH  7 
essentially  only  ammonium  ions  are  present,  at  pH 
12  only  dissolved  ammonia  gas.  This  dissolved  gas 
can  be  released  from  the  liquid  under  proper  con¬ 
ditions,  and  this  is  the  purpose  of  air-stripping  op¬ 
erations  in  mechanical  treatment  plants.  Typical¬ 
ly,  the  pH  is  elevated  with  lime,  the  sludge  is  al¬ 
lowed  to  settle,  and  then  the  clarified  effluent  is 
aerated  or  allowed  to  cascade  over  a  packed  bed  in 
a  stripping  tower.  At  the  high  hydraulic  loading 
rates  commonly  employed  the  concept  is  only  ef¬ 
fective  in  warm  weather  because  of  the  tempera¬ 
ture  constraints  illustrated  in  Figure  1 . 
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Figure  I.  Ammonia-ammonium  ion  distribution 
in  water — influence  of  pH  and  temperature'1. 


2.  Pond  system  hydraulics 

It  is  common  in  the  design  and  analysis  of  sus¬ 
pended  growth  biological  wastewater  treatment 
systems  to  characterize  the  system  as  either  com¬ 
pletely  mixed,  plug  flow,  or  intermediate.  The 
complete-mix  approach  assumes  that  the  reactor  is 
continuously  and  completely  mixed,  so  that  a  sam¬ 
ple  taken  from  any  point  at  a  particular  time 
would  be  identical  to  a  sample  from  any  other 
point  in  the  reactor.  Aerated  systems  with  high-in¬ 
tensity  agitation  approach  these  conditions,  and 
the  equation  describing  biochemical  removal  of  a 
constituent  typically  takes  the  form: 

Af  _  1 

A j  1  +  Kct 

where  A{  =  final  concentration  of  constituent 
/lj  =  initial  concentration  of  constituent 
Kc  =  first  order,  temperature-dependent, 
complete-mix  reaction  rate  constant 
t  =  hydraulic  residence  time. 

Activated  sludge  systems  and  complete-mix  pond 
systems  are  designed  based  on  this  approach. 
Many  partial-mix  (see  Table  1)  ponds  have  also 
been  designed  this  way,  even  though  actual  condi¬ 
tions  do  not  always  fit  the  design  assumptions. 

The  plug  flow  design  assumes  that  a  unit  of 
wastewater  retains  its  identity  throughout  its  pass¬ 
age  through  the  system  and  is  not  mixed  with 
other  increments.  A  sample  taken  near  the  end  of 
the  system  would  not  be  identical  in  character  to  a 
sample  taken  near  the  influent  point.  It  is  there¬ 
fore  not  appropriate  to  compare  input  and  output 
samples  taken  at  the  same  time  and  conclude  that 


a  certain  removal  has  occurred  unless  the  input 
concentration  does  not  vary  with  time.  Since  this 
is  not  the  case  in  typical  systems,  the  true  removal 
can  only  be  determined  by  comparing  the  charac¬ 
teristics  of  the  final  sample  at  discharge  with  its 
own  characteristics  when  it  entered  the  system. 
The  equation  describing  the  biochemical  removal 
of  a  constituent  under  plug  flow  conditions  typi¬ 
cally  takes  the  form: 


where  Kp  =  first-order,  temperature-dependent, 
plug  flow  reaction  rate  constant  and  other  terms 
are  as  defined  previously. 

Pond  systems  that  do  not  have  mechanical  mix¬ 
ing  or  aeration  (i.e.  oxidation  ponds,  facultative 
ponds,  etc.)  are  generally  considered  to  be  plug 
flow  reactors.17  Wind  action  and  thermally  in¬ 
duced  density  currents  can  induce  both  vertical 
and  lateral  mixing,  but  these  are  usually  trans: 
effects  compared  to  the  long  detention  time  i  e 
system.  Bowles  and  Middlebrooks2  have  s  n 
that  short-circuiting  of  flows  occurs  in  many  d 
systems  due  to  configuration,  inlet  and  out 
cation,  etc.  The  actual  detention  time  in  the >n 
studied  in  their  work  ranged  from  29%  to  89%  of 
the  assumed  design  detention  time.  These  reduc¬ 
tions  are  due  to  “dead  spots”  in  the  pond,  so  the 
plug  flow  assumption  would  still  apply  to  the  bulk 
of  the  liquid  moving  through  the  system. 

3.  Pond  biochemistry 

The  ecosystem  in  plug  flow  wastewater  treat¬ 
ment  ponds  is  very  dynamic.  It  is  similar  in  many 
respects  to  the  ecosystem  found  in  natural  ponds, 
but  the  diversity  of  species  tends  to  be  less  and  the 
density  of  those  organisms  present  tends  to  be 
greater  than  the  natural  case  because  of  the  higher 
organic  and  nutrient  loadings  entering  the  system. 
Most  of  the  physical,  chemical  and  biochemical 
reactions  taking  place  in  ponds  have  an  impact  on 
the  nitrogen  present  in  the  system. 

Most  of  the  nitrogen  entering  the  system  will  be 
in  the  ammonia  form,  either  as  a  dissolved  gas  or 
as  ammonium  ions.  The  remainder  will  be  organic 
nitrogen  associated  with  the  solid  particulate  mat¬ 
ter.  Some  of  the  latter  can  settle  directly  to  the 
bottom  due  to  gravity  forces.  Either  ammonia  or 
organic  N  can  be  consumed  by  animals  or  bacteria 
and  then  eventually  settle.  The  ammonia  forms 
can  be  taken  up  by  plants  and  algae  and  also  even¬ 
tually  settle.  These  are  the  major  pathways  sug¬ 
gested  by  Ferrara  and  Avci‘  for  nitrogen  removal 


Table  2.  Data  summary  from  EPA  lagoon  studies.1  “  “  11 


Liquid  temp  (°C) 

Average  algal  cone 
(no./mL) 

Influent  total  N 
(mg/L) 

Nitrogen  removal 
<*) 

Detention 

Ann 

Apr - 

Oct- 

Ann 

Apr- 

Oct- 

Ann 

Apr- 

Oct- 

Ann 

Apr- 

Oct- 

Location 

time,  d* 

avg 

Sep 

Mar 

avg 

Sep 

Mar 

avg 

Sep 

Mar 

avg 

Sep 

Mar 

Peterborough 

107 

n 

18 

5 

398,000 

216,000 

563,000 

17.8 

15.4 

20.1 

43 

50 

38 

Kilmichael 

214 

18.4 

24.2 

13.4 

429,600 

196,000 

580,000 

35.9 

41.8 

30.1 

80 

83 

76 

Eudora 

231 

14.7 

02 

7.2 

412,178 

440,000 

384,000 

50.8 

52.2 

49.3 

82 

80 

83 

Corrine 

42 

10 

14.8 

4.3 

781,000 

1,223,087 

423,072 

14.0 

11.1 

16.9 

46 

56 

36 

'Note:  Peterborough,  Kilmichael  and  Eudora  are  three-cell  systems.  Corrine  is  a  seven-cell  system.  Detention  times,  temperatures,  etc. 
are  given  for  the  first  three  cells  at  Corrine. 


from  pond  systems.  It  is  believed  that  these  can  be 
significant  pathways,  but  there  is  also  a  strong  po¬ 
tential  for  reintroduction  of  nitrogen  into  the 
water  column  from  this  settled  material.  These 
biological  cells  might  contain  up  to  8%  nitrogen, 
but  most  of  the  cellular  material  is  susceptible  to 
decomposition.  It  is  estimated1 1  that  the  nonaero- 
bically-biodegradable  portion  might  range  from 
20%  to  50%.  This  means  that  up  to  80%  of  the 
settled  nitrogen  is  again  potentially  available  for 
solution  in  the  water  column.  Most  of  this  will 
again  be  in  ammonia  form,  and  any  portion  of 
that  in  the  gaseous  form  is  potentially  available 
for  loss  to  the  atmosphere.  The  portion  that  re¬ 
mains  temporarily  with  the  bottom  sediments  is 
not  absolutely  permanent  due  to  the  periodic  re¬ 
suspension  of  some  bottom  material  during  sea¬ 
sonal  overturns,  the  consumption  of  the  suspend¬ 
ed  material  by  protozoa  and  other  animals,  and 
the  anaerobic  decomposition  of  the  material  in- 
place. 

Some  qualitative  indications  regarding  the  role 
of  the  algae  in  the  removal  of  nitrogen  result  from 
an  examination  of  the  data  collected  at  the  four 
facultative  ponds  in  the  recent  EPA  studies.1 10  “II 
Table  2  summarizes  pertinent  average  data  from 
these  four  systems  and  presents  annual  and  sea¬ 
sonal  values.  The  nitrogen  removals  shown  were 
calculated  using  a  plug  flow  assumption  (i.e.  the 
effluent  leaving  the  pond  today,  in  a  30-day  pond, 
is  not  compared  to  the  influent  entering  today,  but 
rather  to  the  influent  that  entered  30  days  ago). 
The  pond  systems  studied  represent  a  wide  geo¬ 
graphical  distribution,  and  all  show  significant 
seasonal  trends  with  respect  to  liquid  temperature 
in  the  system.  However,  the  average  algal  concen¬ 
trations  do  not  show  a  consistent  seasonal  trend, 
although  three  of  the  four  do  show  a  significant 
seasonal  difference  for  nitrogen  removal. 

If  algae  constitute  the  major  factor  responsible 
for  permanent  nitrogen  removal,  it  then  seems 


likely  that  maximum  removals  should  be  observed 
during  or  very  soon  after  peak  algal  growth  peri¬ 
ods.  It  is  possible  to  test  this  hypothesis  by  exam¬ 
ining  the  average  monthly  data  from  these  same 
four  systems.  Figures  2-5  present  the  results  of 
this  analysis.  The  deviation  from  the  mean  annual 
values  is  plotted  for  both  algal  concentration  and 
nitrogen  removal  to  more  clearly  illustrate  the 
trends. 

The  three  pond  systems  in  the  colder  climates 
(Peterborough,  New  Hampshire,  Eudora,  Kansas, 
and  Corinne,  Utah)  exhibit  an  algal  peak  in 
spring,  which  is  probably  due  to  the  seasonal  over¬ 
turn  of  the  pond  and  the  resuspension  of  benthic 
nutrients.  Two  of  these  (Peterborough  and  Co¬ 
rinne)  also  show  an  algal  peak  in  fall  and  early 
winter  when  another  overturn  would  be  expected. 
(The  temperature  must  gradually  change  and  then 
remain  at  about  4°C  for  an  overturn  to  occur. 
These  conditions  did  not  exist  in  the  fall  and  early 
winter  at  the  Eudora  site  during  the  study  period, 
so  only  the  spring  overturn  was  observed.) 

Higher  nitrogen  removals  are  demonstrated  at 
ail  four  sites  during  the  warmer  months  (Eudora  is 
marginal  in  this  respect),  but  there  is  no  apparent 
correlation  between  the  periods  of  higher  nitrogen 
removals  and  the  peak  periods  for  algal  growth. 
This  would  seem  to  indicate  that  algal  uptake  is 
not  the  major  factor  responsible  for  permanent  ni¬ 
trogen  removal. 

The  bacterial  oxidation  of  ammonia  to  nitrate  is 
thought  to  be  minimal  in  pond  systems,  as  indicat¬ 
ed  previously.  The  average  nitrate  concentrations 
in  the  same  four  pond  systems'  ",l6!l  are  summa¬ 
rized  in  Table  3.  Nitrate  concentration  in  the  efflu¬ 
ents  is  very  low  in  all  of  the  systems  and  represents 
less  than  1%  of  the  total  nitrogen  that  entered. 
There  are  no  consistent  seasonal  trends  indicated 
by  the  data,  and  one  system  (Corinne)  shows  a  net 
decrease  in  nitrate  on  a  year-round  basis.  These 
low  values  do  not  prove  that  nitrification  is  not  oc- 
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Figure  2.  Alga I  concentration  and  nitrogen  re¬ 
moval,  Peterborough,  New  Hampshire. 


Figure  3.  Algal  concentration  and  nitrogen  re¬ 
moval,  KUmichael,  Mississippi. 


Figure  4.  Algal  concentration  and  nitrogen  re¬ 
moval,  Eudora,  Kansas. 


Figure  5.  Algal  concentration  and  nitrogen  re¬ 
moval,  Corinne,  Utah. 


curring  in  ponds.  It  may  be  possible  that  signifi¬ 
cant  nitrification  is  taking  place,  but  the  nitrates 
are  then  taken  up  by  the  algae  and/or  denitrified 
and  lost  as  nitrogen  gas.  However,  it  appears  more 
likely  that  nitrification  rates  are  low  in  pond  sys¬ 
tems  since  the  responsible  organisms  require  oxy¬ 
gen  and  relatively  low  levels  of  biodegradable  or¬ 
ganics,  and  are  thought  to  prefer  attachment  to 
surfaces.  These  three  factors  do  not  simultaneous¬ 
ly  exist  to  a  significant  degree  anywhere  in  a  typi¬ 
cal  pond.  (High  levels  of  nitrification  were  ob¬ 
served  in  an  aquaculture  pond  in  Hercules,  Cali¬ 
fornia,  where  suspended  plastic  strips  were  in¬ 
cluded  as  a  support  medium  for  biological 
growth.10) 

If  significant  nitrification  is  occurring  in  typical 
ponds,  it  can  only  be  an  intermediate  ,.age,  as  evi¬ 


denced  by  the  low  nitrate  concentrations  in  the 
system  effluents.  The  permanent  removal  of  ni¬ 
trate  would  then  require  biological  uptake  and  set¬ 
tling  and/or  denitrification.  Denitrification  re¬ 
quires  an  anaerobic  environment  and  an  adequate 
carbon  source.  Both  of  these  factors  will  exist  at 
the  pond  bottom,  but  the  presence  of  significant 
nitrates  at  the  same  location  seems  unlikely.  Par¬ 
ticulate  matter  and  dead  vegetation,  etc.,  may  set¬ 
tle  and  become  part  of  the  benthic  deposit,  but  the 
nitrogen  tied  up  in  those  materials  would  be  in  the 
organic  form,  and  their  subsequent  oxidation  all 
the  way  to  nitrate  seems  unlikely  in  the  anaerobic 
environment  present  on  the  bottom. 

If  algal  uptake/settling,  nitrification,  and  deni¬ 
trification  are  eliminated  as  major  pathways  for 
permanent  nitrogen  removal,  the  only  possibility 
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Table  3.  Nitrates  in  facultative  ponds. 

Nitrite  plus  nitrate  (mg/L) 


Location 

Influent, 

annual 

average 

Effluent * 
Annual  April- 

average  Sept 

Oct- 

Nov 

Peterborough 

0.10 

0.19 

0.22 

0.17 

Kilmichael 

0.11 

0.19 

0.13 

0.25 

Eudora 

<0.1 

0.67 

0.55 

0.78 

Corinne 

0.20 

0.04 

0.04 

0.04 

•This  is  final  effluent  from  the  three-cell  systems  at  Peter¬ 
borough,  Kilmichael  and  Eudora,  and  is  effluent  from  the 
third  cell  at  Corinne. 


left  is  air  stripping  or  volatilization  of  ammonia. 
The  equilibrium  equation  for  ammonia  in  water  is 
represented  by: 

NH/—  NHj  t  +  H  +  . 

Figure  1  illustrates  the  relationship  at  various 
temperatures.  At  typical  ambient  temperatures 
and  low  hydrogen  ion  concentrations  (pH  10.3 
+ ),  essentially  all  of  the  ammonium  is  in  gaseous 
form.  A  high  pH  is  achieved  in  mechanical  treat¬ 
ment  plants  by  adding  lime  or  similar  chemicals. 
Elevated  pH  values  are  also  achieved  in  ponds  due 
to  the  interactions  of  algae  and  the  carbonate  buf¬ 
fering  system.  The  latter  is  represented  by  the  fol¬ 
lowing  equations: 

C02  +  H20  =  H2CO,  =  HCOj"  +  H  + 
HCOr  =CO,  =  +  H  + 

CO,  =  +  H,0  =  HCOr  +  OH‘. 

During  algal  photosynthesis  carbon  dioxide 
(C02)  is  removed  from  solution.  The  result  is  to 
decrease  both  the  bicarbonate  alkalinity  and  the 
hydrogen  ion  concentration  (increases  pH).  Ex¬ 
treme  diurnal  fluctuations  in  the  pH  are  observed 
in  many  ponds  because  of  this  relationship  be¬ 
tween  pH  and  algal  photosynthesis.  On  warm, 
sunny  days  pH  values  exceeding  10-11  are  com¬ 
mon  in  the  surface  layers  of  many  ponds.  Oxygen 
is  another  product  of  algal  photosynthesis  since 
the  carbon  is  extracted  from  the  C02  and  the  ele¬ 
vated  pH  is  typically  accompanied  by  supersatu¬ 
rated  oxygen  conditions.  King1’  has  shown  that 
the  rates  of  solution  of  atmospheric  carbon  diox¬ 


ide  are  too  slow  to  support  optimum  photosynthe¬ 
sis,  so  the  algae  draw  on  the  C02  present  in  the  al¬ 
kalinity  system.  This  suggests  that  bicarbonate  al¬ 
kalinity  has  an  influence  on  the  air  stripping  of 
ammonia  and  that  wastewaters  low  in  alkalinity 
might  not  be  as  responsive  in  pond  systems. 

A  pH  of  least  10-12  is  considered  essential  for 
efficient  Tunonia  air  stripping  in  mechanical  sys¬ 
tems.  The  near-surface  bulk  liquid  in  ponds  may 
only  approach  that  level  for  a  brief  period,  and 
then  only  on  warm,  sunny  days.  However,  the  de¬ 
tention  time  in  the  mechanical  system  is  measured 
in  terms  of  minutes  or  hours,  while  in  the  pond  it 
is  many  days  or  months.  Even  if  the  pH  is  only  8 
or  slightly  above,  and  the  gaseous  fraction  of  am¬ 
monia  stays  less  than  10%,  the  long  detention  time 
can  provide  a  sufficient  period  for  effective  re¬ 
moval.  Even  if  the  pH  of  the  bulk  liquid  is  close  to 
neutral,  it  seems  possible  that  there  could  be 
microsites  at  the  surface  with  an  elevated  pH  due 
to  algal  activity  resulting  in  gaseous  losses  over  the 
long  term. 

The  long  detention  period  also  provides  enough 
time  for  the  nitrogen  entering  a  pond  system  to  cy¬ 
cle  through  a  number  of  oxidation  states  and/or 
cycles  of  plant  uptake  and  decomposition.  Every 
time  a  cycle  is  completed  the  nitrogen  goes 
through  the  ammonia  form  and  a  portion  is  avail¬ 
able  for  loss  as  a  gas.  It  seems  very  unlikely  that 
the  nitrogen  entering  a  pond  will  remain  in  the 
same  discrete  form  throughout  its  passage  through 
the  system.  It  is  far  more  likely  that  one  or  more 
transformations  will  occur  so  that  the  nitrogen 
forms  leaving  a  pond  are  a  result  of  the  activity  in 
the  pond  and  not  a  residue  of  the  incoming  materi¬ 
al.  Therefore,  it  seems  inappropriate  to  compare  a 
particular  form  in  the  effluent  to  the  same  materi¬ 
al  in  the  influent  and  conclude  that  some  absolute 
level  of  removal  has  occurred. 

For  engineering  purposes  it  seems  both  prudent 
and  conservative  to  consider  all  of  the  total  nitro¬ 
gen  entering  and  leaving  the  system,  and  then  to 
assume  that  essentially  all  of  the  nitrogen  leaving 
the  system  is  still  in  the  ammonia  form.  There  is  a 
consensus  among  recent  investigators  that  nitro¬ 
gen  removal  in  ponds  is  dependent  on  at  least  pH, 
temperature  and  detention  time,  even  though 
there  is  dispute  over  the  specific  mechanism  in¬ 
volved.  In  the  opinion  of  this  author,  the  volatili¬ 
zation  pathway  seems  more  likely,  but  there  is  no 
absolute  proof  available.  Fortunately,  it  is  not  ne¬ 
cessary  to  resolve  the  issue  in  order  to  develop  a 
rational  design  procedure. 


DEVELOPMENT  OF  MATHEMATICAL 
MODEL 

1.  General  consideration  and  data  base 

Based  on  the  factors  discussed  above,  a  first 
order  plug  flow  equation  was  adopted  as  the  ap¬ 
propriate  model  to  describe  nitrogen  removal  in 
pond  systems.  The  equation  takes  the  general 
form: 

Ne/N0  =  e-W.t.pH)  (1) 

where  Ne  =  total  nitrogen  in  pond  effluent  at  the 
end  of  the  detention  period  (mg/L) 
N0  =  total  nitrogen  in  pond  influent  at  the 
start  of  the  detention  period  (mg/L) 
K  =  first-order  reaction  rate  constant,  a 
function  of  temperature  (7),  time 
(f),  and  pH. 

The  data  used  to  develop  this  model  were  ob¬ 
tained  from  the  four  EPA  pond  studies  previously 
mentioned:'  10  “J1 

Peterborough,  New  Hampshire.  The  system  is 
a  three-cell  facultative  pond,  operated  in  series, 
with  a  continuous  discharge.  During  the  study  pe¬ 
riod  the  mean  flow  was  101 1  cubic  meters  per  day 
(0.27  million  gallons  per  day).  Based  on  this  mean 
flow  the  hydraulic  detention  period  was  107  days. 
The  design  depth  in  all  cells  was  1.2  m  (4  ft),  and 
the  total  surface  area  for  the  three  cells  was  8.5 
hectares  (21  acres).'1 

Kitmichael,  Mississippi:  A  three-cell  facultative 
pond,  operated  in  series,  with  a  continuous  dis¬ 
charge.  During  the  study  period  the  mean  flow 
was  280  m’/d  (0.07  mgd).  Based  on  this  mean  flow 
the  hydraulic  detention  period  was  214  days.  The 
design  depth  was  about  2  m  (6.6  ft)  in  all  cells,  and 
the  total  surface  area  for  the  three  cells  was  3.3  ha 
(8.1  acres).  '* 

Eudora,  Kansas:  A  three-cell  facultative  pond, 
operated  in  series,  with  a  continuous  discharge. 
During  the  study  period  the  mean  flow  was  500 
m’/d  (0.13  mgd).  Based  on  this  mean  flow  the  hy¬ 
draulic  detention  period  was  231  days.  The  design 
depth  was  1.5  m  (5  ft)  and  the  total  surface  area 
for  the  three  cells  was  7.8  ha  (19.3  acres).10 

Corinne,  Utah:  A  seven-cell  facultative  pond, 
operated  in  series,  with  a  continuous  discharge. 
Only  the  first  three  cells  were  utilized  in  this 
analysis.  Based  on  the  mean  flow  of  690  m!/d 
(0.18  mgd)  during  the  field  study  the  detention 
period  in  the  first  three  cells  was  42  days.  The  de¬ 
sign  depth  was  1.2  m  (4  ft)  and  the  surface  area 
was  2.3  ha  (5.75  acres)  in  these  first  three  cells.1 


All  four  systems  receive  raw  municipal  (essen¬ 
tially  domestic)  wastewater  from  the  adjacent 
community.  Data  collection  continued  for  a  full 
year  at  each  site  to  include  all  seasonal  influences. 
An  ice  cover  was  noted  at  both  Peterborough  and 
Corinne  during  some  winter  months.  The  aver¬ 
aged  data  values  used  for  this  analysis  are  tabu¬ 
lated  in  Appendix  A.  The  nitrogen  concentrations 
given  are  total  nitrogen  values,  and  are  obtained 
by  summing  the  TKN  (total  Kjeldahl  nitrogen) 
and  the  NOj-NOj  (nitrite-nitrate)  data.  Since  the 
TKN  includes  both  the  organic  N  and  the  am¬ 
monia  N  fractions,  it  is  a  close  approximation  of 
the  total  nitrogen  present  in  untreated  wastewater. 

The  effluent  nitrogen  values  listed  in  the  tables 
in  Appendix  A  are  the  actual  average  monthly  val¬ 
ues  measured  for  the  month  indicated.  All  of  the 
other  data  listed  (pH,  temperature,  influent  nitro¬ 
gen,  etc.)  were  obtained  by  applying  the  plug  flow 
assumption  to  the  raw  data.  For  example,  in  the 
42-day  detention  time  pond  at  Corinne,  the  efflu¬ 
ent  nitrogen  listed  for  January  is  the  average  value 
for  the  period  1-31  January.  However,  the  influ¬ 
ent  values  listed  represent  the  characteristics  of  the 
January  effluent  when  it  entered  the  system,  so 
they  are  an  average  of  the  influent  measurements 
for  the  period  20  Nov-20  Dec.  Similarly,  the  tem¬ 
perature  and  pH  values  listed  are  the  average  of 
measurements  made  in  the  pond  during  that  par¬ 
ticular  42-day  detention  period.  For  this  iteration 
no  adjustments  were  made  in  the  detention  times 
used,  to  compensate  for  variations  in  incoming 
flow  and/or  water  losses  due  to  seepage  and  evap¬ 
oration.  These  factors  could  be  very  significant 
where  sustained  seasonal  peak  flows  occur  or  in 
arid  climates  with  a  high  evaporative  loss. 

Each  monthly  value  tabulated  in  Appendix  A 
represents  the  average  of  at  least  five  to  eight 
measurements  made  during  the  period  of  concern. 
Since  the  ponds  all  had  multiple  cells  and  the  in¬ 
puts  and  outputs  were  measured  for  each,  it  is  pos¬ 
sible  to  arrange  the  data  in  a  number  of  configura¬ 
tions  for  analysis.  Table  4  lists  the  cases  selected 
for  development  of  the  model.  The  total  data  base 
for  each  parameter  is  at  least  1000  points,  since 
each  case  has  1 2  monthly  values  and  each  value  is 
an  average  of  five  to  eight  measurements.  The  de¬ 
tention  times  represented  range  from  27  to  231 
days,  the  temperature  from  1.2°  to  26.3  °C,  and 
pH  from  6.4  to  9.5.  It  was  possible  to  identify  spe¬ 
cific  periods  when  the  other  variables  were  rela¬ 
tively  constant  for  analysis  of  the  effect  of  deten¬ 
tion  time,  temperature  and  pH.  With  one  excep¬ 
tion  (Case  9),  all  of  the  data  sets  represent  treat¬ 
ment  performance,  starting  with  raw  untreated 


Table  4.  Data  configurations  used  for  model  development. 


Peterborough 

Case  1  (Table  Al)  AH  three  cells,  t  =  107  days,  inpultocell  I, 
output  form  cell  3,  pH  and  T  are  average  of  all  three  cells. 
Case  2  (Table  A2)  First  cell,  t  -  44  days,  input  to  cell  1,  out¬ 
put  from  cell  1,  pH  and  T  are  average  of  cell  1. 

Case  3  (Table  A3)  First  two  cells,  t  =  74  days,  input  to  cell  I , 
output  from  cell  2,  pH  and  T are  average  of  cells  1  and  2. 


Kilmichael 

Case  4  (Table  A4)  All  three  cells,  /  =  214  days,  input  to  cell  1, 
output  from  cell  3,  pH  and  T  are  average  of  all  three  cells. 
Case  5  (Table  A5)  First  cell,  t  -  136  days,  input  to  cell  1,  out¬ 
put  from  cell  1,  pH  and  T  are  average  of  cell  1. 

Case  6  (Table  A6)  First  two  cells,  T  =  191  days,  input  to  cell 
1,  output  from  cell  2,  pH  and  T  are  average  of  cells  1  and  2. 


sewage  since  that  is  the  purpose  of  the  model. 
Case  9  was  developed  for  comparison  with  Cases  2 
and  1 1  so  there  would  be  additional  data  at  essent¬ 
ially  the  same  detention  time  for  analysis  of  pH  in¬ 
fluences. 

2.  Detention  time  and  temperature 

Selected  data  were  organized  into  groups  that 
reflected  similar  pH  and  temperature  conditions. 
Four  such  groups  could  be  developed  at  4  °C,  8  °C, 
18°C  and  24  °C  respectively.  Figure  6  illustrates 
the  18  °C  group,  which  had  a  pH  range  from  8.3  to 
8.4.  The  equation  of  best  fit  takes  the  form: 

N/Nq  =  c~K7i,  +  A)  (1) 

where  KT  =  temperature-dependent  rate  constant 
t  =  detention  time  (days) 

A  =  constant,  related  to  pH. 

If  nitrogen  removal  were  only  dependent  on  de¬ 
tention  time  then  A  would  be  equal  to  zero  and 
jVj/Aq  would  equal  1  when  t  equals  zero.  The 
slope  of  the  line  in  Figure  6  accurately  describes 
the  time-dependent  rate  (correlation  coefficient 
0.98,  significant  at  0.01  level).  The  time  depend¬ 
ence  is  described  by: 

Nt/N0  =  e‘*W  (2) 

The  rate  constants  at  4  °C,  8  °C  and  24  °C  were 
also  determined  as  described  above,  and  are  listed 
in  Table  5  and  plotted  in  Figure  7. 

The  general  expression  used  to  describe  the  tem¬ 
perature  dependence  of  biological  wastewater 
treatment  processes  is: 


Eudora 

Case  7  (Table  A7)  All  three  cells,  r  =  231  days,  input  to  cell 
1,  output  from  cell  3,  pH  and  T  are  average  of  all  three  cells. 
Case  8  (Table  A8)  First  cell,  /  =  93  days,  input  to  cell  I,  out¬ 
put  from  cell  1 ,  pH  and  T  are  average  of  cell  1 . 

Case  9  (Table  A9)  Second  cell,  f  =  44  days,  input  to  cell  2, 
output  from  cell  2,  pH  and  T  are  average  of  cell  2. 

Case  10  (Table  A10)  First  two  cells,  r  =  137  days,  input  to  cell 
I,  output  from  cell  2,  pH  and  T art  average  of  cells  1  and  2. 

Corinne 

Case  1 1  (Table  Al  1)  All  three  cells,  t  =  42  days,  input  to  cell  1 , 
output  from  cell  3,  pH  and  T  are  average  of  all  three  cells. 
Case  12  (Table  A12)  First  cell,  /  =  27  days,  input  to  cell  1,  out¬ 
put  from  cell  I,  pH  and  T are  average  of  cell  1. 

Case  13  (Table  AI3)  First  two  cells,  t  =  34.5  days,  input  to  cell 
I,  output  from  cell  2,  pH  and  T  are  average  of  cells  1  and  2. 


Kt  =  (3) 

where  KT  =  rate  constant  at  temperature  T 
8  =  temperature  coefficient 
T  =  temperature  ( °C). 

A  regression  analysis  of  the  data  in  Figure  7 
produces: 

Kt  =  0.0064  (1.039)<r-20>  (4) 

where  K}0  =  0.0064 
0  =  1.039 

(correlation  coefficient  0.94,  at  0.06  level). 

The  general  equation  for  time  dependence  be¬ 
comes: 

Nt/N0  =  e-00064f  (at  20  °C).  (5) 

Equation  4  can  then  be  used  to  determine  the 
appropriate  value  of  KT  at  temperatures  other 
than  20  °C. 


Table  5.  /ft(T)  constant 
vs  temperature. 

Temperature  f°C) _ K_ 

4  0.0031 

8  0.0050 

18  0.0065 

24  0.0068 
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Figure  6.  Ne/N0  vs  detention  time.  Data 
from  Kilmichael,  Eudora  and  Corinne,  for 
similar  temperature  (18°C)  and  pH  condi¬ 
tions  (8. 3-8.4). 


9.  pH  dependence 

It  is  possible  to  determine  the  pH  influence  in¬ 
dependently  by  comparing  data  groups  at  essen¬ 
tially  the  same  detention  time  and  temperature 
(Cases  2,  9  and  11).  Figure  8  illustrates  a  typical 
set  at  an  average  temperature  of  16.5  °C.  These 
data  can  be  described  by: 


N'/N0  =  e-0.313(pH-6.6>. 


(6) 


The  correlation  coefficient  was  0.86  and  is  sig¬ 
nificant  at  the  0.05  level,  indicating  that  a  good 
approximation  is  provided. 

A  temperature  dependence  was  indicated  by  the 
data,  so  rate  constants  were  determined  for  other 
temperature  conditions  within  this  42-44  day  de¬ 
tention  time  group.  These  results  are  plotted  in 
Figure  9,  and  the  data  can  be  described  reasonably 
well  with  the  same  theta  value  used  in  equation  4: 


*PH,„  *  *W0B<*)<M» 


(7) 


where  X_H  =  pH  rate  constant  at  temperature 
m  r(oC) 

XpH  =  pH  rate  constant  at  20 °C 
(20)  =  0.388 
0  =  1.039. 


Therefore,  at  20  °C  eq  6  becomes: 


Figure  7.  Time-dependent  rate  constant  vs  tern-  Figure  8.  pH  vs  N^/Np.  Data  points  from 

perature.  Eudora,  Peterborough  and  Corinne,  t  =  44 

days,  average  temperature  16.5  °C. 


10  20 
Temperature  (*C) 


Figure  9.  pH  rate  constant  vs  temperature. 


This  equation  will  predict  some  nitrogen  remov¬ 
al  down  to  a  pH  of  6.6,  which  is  compatible  with 
the  actual  ammonia  removal  data  from  all  of  the 
sites  investigated.  Nitrogen  removal  at  such  a  low 
pH  must  mean  that  either  some  mechanism  other 
than  volatilization  is  responsible  or  that  there  are 
microsites  at  the  water  surface  with  suitable  pH 
even  when  the  bulk  liquid  gets  down  to  a  pH  of 
6.6. 

4.  Design  equation 

It  is  possible  to  combine  eq  1  and  8,  since  the 
temperature  dependence  is  essentially  the  same  for 
both  rate  constants: 

Ne/N0  =  e-*r<'+-4> 

_  e-[0.0064</)  + 0.388  (pH-6. 6)1 
_  e-0.0064[/  +  60.6(pH-6.6)l 
_  e-  Kt  [f  +  60.6  (pH-6.6)] 

where  KT  =  temperature-dependent  rate  constant 

KT  =  ^20°c  (®) 

Kjqoq  =  0.0064 

0  =  1.039 

and  all  other  terms  are  as  defined  previously. 
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Figure  10.  Predicted  vs  measured  effluent  nitro¬ 
gen  for  all  data  used  in  model  development. 
Data  range:  pH  6. 4-9.5,  temperature  1.2°  to 
26.3  °C,  detention  time  27  to  23!  days. 
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Figure  11.  Nitrogen  removal  in  all  three  cells  at 
Peterborough. 

This  design  equation  was  applied  to  all  of  the 
156  data  sets  and  is  compared  to  the  values  actual¬ 
ly  measured  in  the  field  in  Figure  10.  A  reasonably 
good  fit  and  a  consistent  trend  are  demonstrated. 
Figures  11-13  show  the  predicted  versus  measured 
effluent  values  on  a  monthly  basis  for  the  three 
cases  developed  for  Peterborough,  New  Hamp¬ 
shire.  The  influent  nitrogen  (determined  using 
plug  flow  assumption)  is  shown  in  comparison,  as 
is  the  period  of  continuous  ice  cover. 


Figure  12.  Nitrogen  removal  in  first  two  cells  at 
Peterborough. 


5.  Model  validation 

Not  all  of  the  1S6  data  sets  available  were  used 
in  the  model  development,  so  to  some  degree 
Figure  10  represents  at  least  a  partial  validation  of 
the  design  equation.  However,  it  was  considered 
essential  to  obtain  field  data  from  totally  different 
sites  and  preferably  from  totally  different  geo¬ 
graphical  settings  and  operating  conditions.  This 
is  a  difficult  task  since  the  four  systems'  '• 16 11  used 
represent  the  most  extensive  and  most  complete 
data  available  anywhere  in  the  literature.  After 
considerable  searching,  a  limited  amount  of  data 
from  five  treatment  ponds  were  obtained.  These 
were: 

Gulfport,  Mississippi:  single  cell,  continuous 
discharge,  5-day  detention  time  polishing  pond 
following  aerated  ponds,  design  flow  2900  mJ/d.5 

Drayton  Valley,  Alberta,  Canada:  two  “long” 
detention  cells,  69  days  detention,  5.5  ha,  1.52  m 
deep,  continuous  discharge,  1088  m’/d  flow,  3550 
people.’ 

Stettler,  Alberta,  Canada:  first  “long”  deten¬ 
tion  cell,  75  days  detention  time,  5.34  ha,  1.52  m 
deep,  continuous  discharge,  1082  m’/d  flow,  3800 
people.’ 

Inisfail,  Alberta,  Canada:  first  “long”  deten¬ 
tion  cell,  22.6  days  detention  time,  1.18  ha,  1.52  m 
deep,  continuous  discharge,  800  m’/d  flow,  2600 
people’. 

Lacombe,  Alberta,  Canada:  single  cell,  drained 
in  October,  refill  and  storage  October  to  March, 
discharge  continuous  and  equal  to  input  March  to 


Figure  13.  Nitrogen  removal  in  first  cell  at  Peter¬ 
borough. 


September.  Average  detention  time  during  dis¬ 
charge  period  165  days,  6.73  ha,  1.52  m  deep,  623 
mVd  average  flow,  3140  people.’ 

Pertinent  data  from  these  five  sites  are  tabu¬ 
lated  in  Appendix  B.  The  effluent  values  listed  are 
the  actual  measured  results  for  the  time  indicated. 
The  other  entries  were  determined  using  the  plug 
flow  assumption  previously  described.  The  data 
from  the  four  Canadian  sites  represent  individual 
grab  samples  obtained  during  the  period  1961  to 
1964.  The  values  listed  for  a  particular  month  then 
represent  the  average  of  all  samples  collected  in 
that  same  month  during  the  study  period  (for  ex¬ 
ample,  the  nine  data  entries  in  Table  B3  represent 
18  separate  measurements).  There  was  no  data 
collection  in  some  months  and  they  are  not  listed. 

These  five  sources  yielded  a  total  of  37  data  sets 
that  could  be  used  for  analysis  and  model  verifica¬ 
tion.  The  predicted  nitrogen  was  calculated  using 
eq  9  and  the  data  in  Tables  B1-B5.  These  predicted 
values  are  compared  to  the  actual  effluent  nitro¬ 
gen  concentrations  in  Figure  14.  Both  the  fit  and 
the  trend  are  reasonably  close,  so  it  can  be  con¬ 
cluded  that  the  design  equation  can  be  used  to  give 
a  valid  estimate  of  nitrogen  removal  in  pond  sys¬ 
tems. 

The  design  equation  should  therefore  be  applic¬ 
able  for  any  typical  pond  system  located  in  the 
arctic,  subarctic  or  temperate  zone.  Figure  15 
compares  the  predicted  nitrogen  to  the  actual 
measured  values  for  the  short  detention  pond  in 
Gulfport. 
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Figure  14.  Predicted  vs  measured  effluent  nitro¬ 
gen  for  all  sites  used  for  model  verification. 
Data  range:  pH  6. 9-9. 2,  temperature  1°  to 
28  °C,  detention  time  5  to  165  days. 


Figure  15.  Predicted  vs  measured  effluent  ni¬ 
trogen,  polishing  pond,  Gulfport,  Mississippi 
(X  =5  days,  pH  =  6.9  to  7.6,  T  =  10.9°  to 
28  °C). 
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EFFECT  OF  AN  ICE  COVER 

A  continuous  ice  cover  existed  at  Peterborough 
for  at  least  three  months  and  at  Corinne  for  a  less¬ 
er  period.  However,  the  design  equation  will  still 
predict  a  significant  amount  of  nitrogen  removal 
during  these  same  periods  and  this  is  verified  by 
measured  values  at  both  sites.  Even  if  all  other 
conditions  in  the  liquid  would  allow  the  presence 
of  dissolved  ammonia  gas  it  is  difficult  to  see  how 
the  gas  could  escape  to  the  atmosphere  from  under 
a  continuous  ice  cover.  Some  other  mechanism 
may  be  temporarily  responsible  for  nitrogen  re¬ 
moval  during  periods  of  ice  cover,  with  volatiliza¬ 
tion  again  a  factor  during  ice-free  periods. 

As  might  be  expected,  the  ice  was  quite  thick  at 
the  Canadian  sites  examined  and  approached  1  m 
or  more  in  depth  by  the  end  of  the  winter  (Tables 
B2-B5).  This  is  at  least  Vs  of  the  total  design  depth 
of  the  ponds.  If  about  0.3  m  is  allowed  for  sludge 
accumulation,  this  means  that  less  than  15 Vo  of 
the  design  volume  was  available  to  actually  treat 
the  incoming  wastewater.  This  reduced  detention 
time  coupled  with  the  low  temperatures  is  one  rea¬ 
son  why  ice-covered  ponds  in  northern  climates 
perform  poorly  during  the  extreme  winter  months. 
Effluent  quality  can  be  expected  to  be  little  better 
than  equivalent  to  primary  under  these  conditions. 

Another  water  quality  effect  can  be  observed 
due  to  this  ice  formation  since  most  of  the  dis¬ 


solved  constituents  in  the  wastewater  will  not  be 
incorporated  in  the  ice  crystals  during  the  freezing 
process  but  will  accumulate  in  the  remaining  li¬ 
quid.  This  has  been  observed  during  the  freezing 
of  sea  ice  and  has  been  used  as  a  treatment  process 
for  salt  reduction  in  brackish  waters. 

At  the  Peterborough  site,  the  alkalinity  in  the 
untreated  wastewater  averaged  about  104  mg/L 
on  an  annual  basis  and  did  not  vary  much  from 
month  to  month.  During  the  warm  summer 
months  when  algal  activity  was  greatest  the  alka¬ 
linity  in  the  pond  was  about  84  mg/L.  The  reduc¬ 
tion  is  due  to  the  interaction  of  algae,  CO>  and  the 
bicarbonate  system  previously  discussed.  During 
the  winter  months  the  alkalinity  increased  and  the 
investigators  attributed  the  increase  to  CO,  pro¬ 
duction  by  anaerobic  decomposition.1  However, 
in  February  the  average  alkalinity  in  the  three  cells 
was  164  mg/L,  which  is  about  a  60 Vo  increase  over 
background,  and  the  low  liquid  temperatures 
(3°C)  prevailing  during  the  period  would  not 
favor  significant  anaerobic  decomposition.  An  al¬ 
ternate  or  companion  possibility  might  be  concen¬ 
tration  effects  due  to  rejection  during  ice  forma¬ 
tion. 

Similar,  but  even  more  dramatic  responses  can 
be  observed  at  the  Canadian  systems  used  for 
model  validation  (Tables  B1-B5).  At  Lacombe, 
Alberta,  for  example  (Table  B5),  an  ice  depth  of 
about  0.5  m  is  recorded  for  January  and  an  alka- 
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linity  of  1775  mg/L.  This  is  an  increase  of  about 
325%  over  the  average  value  measured  during  ice- 
free  periods  (575  mg/L).  The  nitrogen  in  this  pond 
also  increased  to  83  mg/L  during  January  as  com¬ 
pared  to  an  average  of  51  mg/L  in  the  incoming 
wastewater  or  an  average  in  the  pond  of  about  16 
mg/L  during  ice-free  periods.  Since  the  0.5  m  of 
ice  represents  about  30%  of  the  pond  volume  it 
seems  likely  that  at  least  some  of  these  increases 
are  due  to  the  concentrating  effects  of  ice  forma¬ 
tion. 

Fertuck*  observed  that  the  ice,  after  natural 
freezing  of  brackish  ponds  in  western  Canada, 
contained  about  20%  of  the  salts  that  were  in  the 
original  water  and  he  developed  a  desalination 
process  to  take  advantage  of  the  phenomenon. 
Adopting  his  20%  factor  as  a  “rule  of  thumb” 
and  further  assuming  that  the  ice  would  have  a 
density  of  about  0.8  g/cm\  it  is  possible  to  calcu¬ 
late  the  concentration  in  the  remaining  liquid  for  a 
particular  thickness  of  ice,  depth  of  pond,  and 
concentration  of  the  same  material  in  the  original 
unfrozen  water. 

It  is  possible  to  avoid  repetitive  calculations  and 
to  develop  a  general  expression  by  normalizing 
both  sides  of  the  equation: 

Concentration  Under  Ice/Original  Unfrozen 
Concentration  =  /  (Ice  Thickness/Pond  Depth) 

or  C  =  /(/) 

where  C  =  concentration  increase  (%) 

/  =  ratio  of  ice  thickness  to  pond  depth. 

When  C  equals  one,  /  will  equal  zero;  when  /  ap¬ 
proaches  1 .0,  C  will  approach  infinity  as  the  con¬ 
centration  increases  in  the  last  drop  of  water  just 
before  it  freezes.  In  fact,  that  last  drop  would 
probably  not  freeze  at  all  due  to  the  increased  salt 
concentration.  It  is  unlikely  even  in  the  most  ex¬ 
treme  climate  that  the  I  factor  would  exceed  0,7  to 
0.8  for  treatment  ponds  with  a  continuous  input 
of  relatively  warm  sewage.  An  equation  can  be  de¬ 
rived  for  the  boundary  conditions  described 
above,  with  an  upper,  practical  limit  for  /  at  0.8: 

C  =  1  -t-  1.72  /*  56  (10) 

where  C  =  concentration  increase  (as  a  decimal) 
=  concentration  in  remaining  liquid 
(mg/L)/original  concentration  in  un¬ 
frozen  liquid  (mg/L) 

I  =  ratio  of  ice  thickness  (as  a  decimal) 

=  ice  thickness/pond  depth. 


0  0.2  0.4  0.6  0.8 

_  /  ice  thickness  \ 

1  \  pond  depth  / 


Figure  16.  Concentration  increase  in  re¬ 
maining  liquid  vs  ice  cover  on  a  pond. 
Data  points  are  from  Lacombe,  Stettler 
and  Drayton  Valley. 


This  expression  is  plotted  in  Figure  16  and  com¬ 
pared  to  the  observed  concentration  increases  for 
alkalinity  at  the  Canadian  sites  included  in  this 
analysis.  The  fit  is  not  statistically  significant  but 
eq  10  can  provide  at  least  a  first  approximation  of 
the  increase  in  concentration  that  might  be  ex¬ 
pected  under  an  ice  cover  in  wastewater  ponds.  It 
should  be  valid  for  any  dissolved  constituent  that 
remains  dissolved  (i.e.  no  settling,  precipitation, 
etc.). 


APPLICATION  OF  THE 
DESIGN  EQUATION 

Input  data  required  for  solution  of  the  design 
equation  include  the  influent  nitrogen  concentra¬ 
tion  (N0),  the  detention  time  in  the  system  (/),  the 
temperature  of  the  liquid  (T),  and  the  pH  of  the  li¬ 
quid  (pH).  Sources  and  methods  for  determining 
these  values  are  described  below.  Since  most  regu¬ 
latory  agencies  specify  30-day  averages  as  limits  it 
is  desirable  to  repeat  the  calculation  to  predict  a 
monthly  average  effluent  nitrogen  as  well  as  the 
annual  average.  The  basic  calculation  should  be 
concerned  with  the  input  to  the  first  cell  and  the 
output  from  the  final  cell.  Three  cells  are  typically 
used  in  most  current  designs  and  the  minimum  de¬ 
tention  time  (0  is  that  necessary  to  satisfy  the 
BOD  requirements  and  is  determined  separately. 

1.  Determination  of  input  nitrogen 

Actual  data  should  be  used  whenever  available, 
or  data  from  similar  facilities.  As  shown  by  the  ac- 


tual  data  compiled  in  Appendices  A  and  B  the  in¬ 
put  nitrogen  concentration  varies  significantly 
over  a  year’s  time  in  typical  municipalities.  If  data 
are  not  available,  the  assumption  of  a  constant 
average  value  will  be  necessary.  As  indicated  pre¬ 
viously,  the  nitrogen  in  raw  domestic  sewage 
ranges  to  over  50  mg/L.  The  following  values  are 
suggested  as  preliminary  estimates: 

Weak  Medium  Strong 

sewage  sewage  sewage 

_ (BOD  -120)  (BOD  «  220)  (BOD  -420) 

Total  N:  20  35  60 

If  actual  monthly  averages  are  available  then  the 
input  values  for  calculation  should  be  determined 
using  the  plug  flow  approach  previously  de¬ 
scribed. 

2.  Determination  of  pond  temperature 

The  temperature  in  a  wastewater  pond  is  a  func¬ 
tion  of  pond  geometry,  detention  time,  influent 
temperature  and  air  temperature.  The  major  heat 
losses  are  from  the  water  surface  by  convection 
and  radiation.  Eckenfelder4  developed  the  follow¬ 
ing  equation  for  estimating  the  bulk  liquid  temper¬ 
ature  of  wastewater  treatment  ponds: 

Tw  =  (T0  +  kTA)/(l+k)  (11) 

where  Tw  =  mean  temperature  in  pond  ( °F) 

7)  =  influent  temperature  ( °F) 

7"a  =  air  temperature  ( °F) 
k  =  proportional  factor 
=  1.6  t/d 
t  =  time  (days) 
e  =  pond  depth  (ft). 

It  is  necessary  to  use  U.S.  customary  units  in 
solving  the  equation  since  the  proportionality  fac¬ 
tor  includes  a  number  of  coefficients  and  con¬ 
stants  that  cannot  be  easily  converted.  The  equa¬ 
tion  is  not  valid  after  ice  formation  starts 
(Tv  =  32  °F).  The  equation  was  originally  devel¬ 
oped  for  a  complete-mix  situation  but  provides  a 
reasonably  accurate  estimate  for  plug  flow  condi¬ 
tions  as  well.  Figure  17  compares  the  predicted 
pond  temperature  using  eq  1 1  to  actual  values  for 
the  Gulfport,  Mississippi,  pond.  The  influent  tem¬ 
peratures  (Tj)  can  either  be  obtained  from  records 
or  estimated  by  the  engineer.  The  air  temperatures 
(TA)  can  be  obtained  from  local  weather  records. 
The  first  iteration  should  calculate  the  monthly 
average  water  temperature  during  the  plug  flow 
detention  period  of  concern.  The  average  water 


Figure  17.  Predicted  vs  measured  pond  temper¬ 
ature,  Gulfport. 

temperature  during  the  plug  flow  detention  period 
of  concern  can  then  be  determined  from  this  list. 

3.  Determination  of  pH  in  the  pond 

The  pH  of  pond  contents  will  vary  on  a  daily 
basis  during  warm,  sunlit  periods,  and  will  seldom 
be  equal  to  the  pH  of  the  entering  wastewater  due 
to  the  interaction  between  the  algae  and  the  C02/ 
alkalinity  system  in  the  pond.  If  possible,  actual 
pH  data  should  be  obtained  from  an  operating 
pond  that  is  subject  to  conditions  similar  to  those 
expected  for  the  system  under  design.  The  data  in 
Appendices  A  and  B  can  be  used  for  that  purpose. 
Examination  of  these  data  demonstrates  the  inter¬ 
relationships  of  pH  and  alkalinity.  At  Peter¬ 
borough,  for  example,  where  the  alkalinity  was 
typically  less  than  100  mg/L,  the  pH  remained 
about  7.  In  contrast,  at  Corinne,  where  the  alka¬ 
linity  was  always  above  500  mg/L,  the  pH  was 
typically  9  or  higher. 

The  relationship  between  average  pH  in  the 
overall  system  and  alkalinity  in  the  wastewater 
was  determined  for  the  four  systems  used  in  devel¬ 
opment  of  the  nitrogen  removal  model: 

pH  =  7.3  e  0  0005  Alk  (12) 

where  pH  =  the  average  pH  in  the  pond 
Alk  =  alkalinity  of  the  entering  wastewa¬ 
ter  (mg/L). 

This  equation  is  shown  graphically  in  Figure  18. 
It  will  tend  to  provide  a  conservative  estimate  of 
pH  in  the  pond. 

The  pond  designer  will  usually  initially  deter¬ 
mine  the  minimum  detention  time  required  for 
BOD  removal.  Equation  9  can  then  be  used  to  esti¬ 
mate  the  nitrogen  removal  that  will  occur  under 
the  pH  and  temperature  conditions  (determined  as 
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Figure  18.  pH  in  wastewater  ponds  vs  alkalinity. 
Data  points  are  average  observed  values  from 
Peterborough,  Eudora,  Corinne  and  Kilmichael. 


described  above)  prevailing.  If  further  nitrogen  re¬ 
moval  is  then  specifically  required,  the  designer 
need  only  consider  the  nitrogen  remaining  in  the 
pond  effluent.  Equation  9  can  also  be  used  to  de¬ 
termine  if  additional  detention  time  in  the  pond 
can  satisfy  all  of  the  nitrogen  removal  require¬ 
ments.  The  costs  for  providing  this  additional  time 
can  then  be  compared  to  the  costs  for  alternative 
nitrogen  removal  methods.  For  this  purpose,  eq  9 
takes  the  following  form: 

t  =  In  (A'e/Ay/(-A'T)-60.6  (pH-6.6).  (13) 

All  terms  were  defined  previously. 

Example:  At  a  pH  value  of  9  and  temperature 
of  23  °C  find  the  detention  time  required  for  80% 
nitrogen  removal. 

N^/Nq  =  0.2 

Kt  =  KM6)(T2a)  =  0.0064(1.039)<23-20) 

=  0.0072 

/  =  In  (0.2)/(-0.0072)-60.6  (2.4) 

=  224-145 

=  78  days  for  80%  nitrogen  removal. 

The  design  equation  can  also  be  very  useful  in 
the  design  and  operation  of  land  treatment  sys¬ 
tems,  since  nitrogen  is  usually  the  limiting  factor 
determining  the  land  area  required  for  treatment. 
Under  ideal  conditions  even  a  few  weeks  detention 
in  a  pond  can  have  a  significant  effect  on  nitrogen, 
and  can  therefore  very  significantly  influence  the 
size  and  the  costs  of  the  final  land  treatment  com¬ 
ponent.  Recognition  of  the  nitrogen  changes  de¬ 
scribed  by  this  equation  may  also  allow  selection 
of  more  desirable  crops  or  changes  in  operational 
procedure.  Corn,  for  example,  is  a  high  market 
value  crop,  but  it  has  a  variable  need  for  nitrogen 
during  the  growing  season  (higher  need  in  early 
stages).  It  should  be  possible  to  design  a  system  to 
bypass  part  of  a  pond  during  the  early  part  of  the 
growing  season  and  then  use  final  pond  effluent 


during  the  latter  stages.  This  would  more  closely 
match  the  nitrogen  needs  of  the  crop  and  ensure 
groundwater  protection  at  all  times. 

A  combination  treatment/storage  pond  is  often 
found  at  land  treatment  systems  in  the  northern 
temperate  zone.  Typically  the  treatment  cell  may 
have  1  to  3  days  detention  with  some  aeration  (in 
use,  or  standby)  for  odor  control.  The  detention 
time  in  the  storage  cell(s)  is  determined  using  pro¬ 
cedures  in  the  Corps  of  Engineers  Land  Treatment 
Design  Manual.  “  These  storage  cells  are  usually 
drawn  down  by  late  fall  and  allowed  to  refill  dur¬ 
ing  the  cold  winter  months.  A  similar  concept, 
termed  “controlled  discharge,”  is  commonly  used 
in  Michigan  and  adjacent  states  in  the  north  cen¬ 
tral  United  States.  In  this  approach,  wastewater  is 
held  in  the  pond  during  low  temperature  periods 
and  often  also  during  high  algal  growth  periods, 
and  is  discharged  once  or  twice  per  year  when  con¬ 
ditions  are  appropriate  in  the  receiving  stream. 
The  Lacombe  system  used  for  model  verification 
is  a  pond  of  this  basic  type.’ 

The  geometry  and  hydraulics  of  these  storage 
cells  do  not  completely  satisfy  plug  flow  condi¬ 
tions  during  the  filling  and  storage  periods.  The 
first  increment  of  wastewater  during  the  refilling  is 
likely  to  spread  over  the  entire  pond  bottom.  If 
each  subsequent  increment  is  assumed  to  do  like¬ 
wise  during  the  filling  period,  the  contents  will  ap¬ 
proach  the  complete-mix  case  since  a  vertical  sam¬ 
ple  of  liquid  at  any  point  in  the  pond  will  be  identi¬ 
cal  to  any  other.  It  is  only  when  the  pond  starts  to 
discharge  that  it  can  be  assumed  that  the  incoming 
wastewater  displaces  an  equivalent  volume  at  the 
outlet  and  plug  flow  is  restored. 

As  a  result,  the  calculation  of  nitrogen  loss  from 
these  storage  ponds  is  a  two-step  procedure.  The 
“average”  detention  time  for  use  in  eq  9  to  calcu¬ 
late  nitrogen  losses  during  the  filling  and  storage 
period  is  half  of  the  actual  detention  period.  The 
pH  and  temperature  conditions  are  the  averages 
over  the  entire  period.  This  will  predict  the  nitro- 
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gen  concentration  in  the  pond  just  prior  to  dis¬ 
charge.  Once  discharge  commences,  the  plug  flow 
assumption  prevails.  The  t  is  then  the  total  deten¬ 
tion  time  and  the  nitrogen  previously  calculated  is 
the  N0  for  eq  9  until  the  pond  is  drawn  down  or 
the  contents  displaced  by  fresh  wastewater.  This 
approach  will  allow  a  conservative  estimate  of  ni¬ 
trogen  removal  from  any  combination  of  treat¬ 
ment/storage  with  and  without  bypasses. 
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APPENDIX  A:  POND  DATA  USED  FOR  MODEL  DEVELOPMENT 


Table  A1.  Peterborough,  three  cel  Is,  107  days. 


Final  e? fluent 
total  N  (mg A) 

No 

Input,  total 

N  <mg/U 

Temp 

CO 

pH 

Alkalinity 

(CaC03) 

( mgA) 

Jan 

24.2 

33.6 

4.3 

7.1 

71 

Feb 

30.2 

33.9 

2.7 

7.0 

76 

ffer 

23.1 

2.3 

7.0 

88 

Apr 

20.8 

33.5 

3.0 

7.0 

100 

May 

22.5 

33.5 

8.3 

7.0 

101 

Jun 

16.6 

29.2 

15.0 

7.1 

91 

Jul 

16.7 

27.4 

21.3 

7.1 

85 

Aug 

8.5 

27.6 

23.0 

7.1 

92 

Sep 

7.5 

28.9 

21.3 

7., 

90 

Oct 

11.3 

28.9 

16.6 

7.2 

80 

Nov 

11.9 

28.4 

11.3 

7.2 

72 

Dec 

20.1 

30.3 

6.7 

7.2 

71 

Table  A2.  Peterborough,  first  cel  I,  44  days. 


Np 

Final  effluent 
total  N  (mq/L) 

No 

Input,  total 

N  (mq/L) 

Temp 

CC) 

pH 

Jan 

32.6 

34.5 

4.2 

7.1 

Feb 

32.6 

33.7 

4.9 

6.8 

Mar 

25.5 

34.4 

3.5 

6.8 

Apr 

23.3 

31.7 

3.3 

6.8 

May 

29.0 

31.5 

4.9 

6.5 

Jun 

24.9 

28.6 

11.0 

6.6 

Jul 

25.3 

31.6 

16.8 

6.8 

Aug 

19.2 

28.3 

23.6 

7.0 

Sep 

17.2 

27.3 

22.4 

7.3 

Oct 

31.2 

17.6 

7.3 

Nov 

27.9 

34.2 

14.5 

7.1 

Dec 

27.8 

36.2 

8.7 

7.2 

Table  A3, 
74  days. 

Final 

total 

Peterborough,  first  and  second  cells, 

Np  Nc 

effluent  Input,  total  Temp 

N  (mq/L)  N  (mg/L)  CC)  pH 

Jan 

30.9 

36.6 

4.3 

7.1 

Feb 

33.5 

31.3 

3.6 

6.9 

Mar 

26.1 

34.6 

3.1 

6.8 

Apr 

20.7 

31.6 

3.3 

6.8 

May 

24.3 

31.7 

5.8 

6.5 

Jun 

19.3 

24.9 

15.9 

6.4 

Jul 

17.1 

23.7 

23.6 

7.2 

Aug 

18.5 

31.6 

23.0 

7.0 

Sep 

10.6 

28.9 

21.4 

7.2 

Oct 

18.2 

26.9 

17.2 

7.3 

Nov 

21.2 

31.3 

14.5 

7.2 

Dec 

21.6 

34.2 

5.5 

7.3 

r 

• 

0- 

1 

y 

► 

Table  A4.  Kllmlchael,  three  cel  Is 

,  214  days. 

r 

Ne  % 

Final  effluent  Input,  total  Temp 

Alka  1 1  n  Ity 
(CaC03> 

total  N  (mqA)  N 

(mq/L) 

CC) 

P_H 

( mq/L) 

Jan 

9.1 

41.2 

19.3 

8.3 

131 

Feb 

7.8 

49.1 

17.3 

8.2 

131 

far 

7.5 

47.9 

15.0 

8.1 

127 

Apr 

5.9 

53.5 

13.6 

8.0 

120 

May 

4.7 

40.8 

13.6 

7.9 

112 

Jun 

4.4 

32.9 

14.8 

8.0 

103 

Jul 

8.1 

30.9 

17.1 

8.1 

97 

Aug 

5.2 

24.4 

21.9 

8.3 

98 

Sep 

5.3 

27.3 

21.9 

8.4 

103 

Oct 

3.9 

26.1 

23.1 

8.4 

1 19 

_  Nov 

2.9 

34.9 

23.4 

8.4 

1  19 

Dec 

9.0 

25.1 

22.2 

8.4 

127 

Table  A5.  Kllmlchael, 

f  trst  cel  1 

,  1 36  day  s. 

V 

Final  effluent 

No 

Input,  total 

Temp 

total  N  (mq/_) 

N 

(mqrt.) 

CC) 

PH  . 

Jan  12.0 

50.7 

15.3 

8.0 

Feb  11.6 

47.1 

13.1 

8.0 

Mjr  12.4 

37.1 

11.1 

8.0 

L 

Apr  9.0 

32.1 

11.2 

8.0 

m 

May  9.2 

27.7 

13.7 

8.1 

P 

Jun  7.8 

25.9 

21.0 

8.1 

Jul  8.8 

26.7 

21.0 

8.3 

Aug  9.0 

30.5 

24.7 

8.4 

Sep  9.8 

30.0 

26.3 

8.5 

Oct  7.7 

45.2 

23.5 

8.5 

Nov  14.8 

50.7 

19.7 

8.2 

L_ 

Dec  14.8 

50.7 

19.7 

8.2 

m 

Table  A6.  Kllmlchael,  first  and  second  cells. 

191  days. 

• 

Final  effluent 

No 

Input,  total 

T  emp 

total  N  (mq/L) 

N 

(mq/L) 

CC) 

.  PH 

Jan  10.6 

49.1 

18.6 

8.4 

Feb  8.9 

47.9 

16.5 

8.2 

Mar  8.4 

53.5 

13.4 

8.1 

Apr  6.7 

40.8 

12.5 

8.2 

May  5.7 

33.3 

13.0 

8.2 

Jun  4.5 

30.9 

15.3 

8.3 

Jul  6.6 

24.4 

18.4 

8.5 

Aug  5.8 

27.3 

21.5 

8.6 

Sep  8.0 

26.1 

23.7 

8.8 

Oct  4.6 

34.9 

24.6 

8.7 

Nov  10.4 

41.2 

24.0 

8.7 

Dec  9. 7 

41.2 

21.7 

8.6 

20 
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Table  A7.  Eudora,  three  cells. 

231  days. 

- 

FI  nal 

Ne  No 

effluent  Input,  total  Temp 

Alka 1 1 nlty 
(CBCO3) 

I 

tota  1 

N  (mqA.)  N 

(mq/U  CC) 

pH 

( mg  A.) 

Jan 

9.1 

39.0  17.0 

8.4 

270 

Feb 

9.6 

48.0  13.9 

8.3 

280 

Msr 

9.8 

60.5  10.9 

8.3 

290 

Apr 

12.7 

53.0  8.9 

8.3 

300 

May 

8.7 

45.0  9.7 

8.3 

300 

P 

Jun 

7.7 

47.5  10.0 

8.3 

320 

Ju! 

11.6 

49.0  12.6 

8.4 

300 

Aug 

10.3 

48.5  15.8 

8.5 

280 

Sep 

9.0 

52.0  18.7 

8.6 

280 

Oct 

9.5 

52.0  20.7 

8.6 

270 

Nov 

10.0 

59.0  21.2 

8.5 

260 

Dec 

7.7 

55.5  19.6 

8.5 

260 

i 

Table  A8.  Eudora,  first  cell. 

93  days. 

Final  effluent 

No 

Input,  total 

Temp 

-  «i 

1  , 

total  N  (mq/t) 

N  (mq/L) 

CC) 

PH 

Jan 

23.7 

44.0 

8.0 

7.8 

Feb 

25.7 

51.0 

4.7 

7.8 

Mjt 

21.6 

47.0 

4.2 

7.9 

Apr 

31.1 

50.0 

6.4 

8.1 

May 

29.5 

54.0 

11.1 

8.1 

i 

< 

Jun 

12.9 

50.0 

14.7 

8.2 

Jul 

23.5 

68.0 

21.6 

8.2 

Aug 

25.4 

43.0 

24.3 

8.3 

Sep 

14.0 

35.0 

24.6 

8.2 

Oct 

21.5 

61.0 

22.8 

8.2 

Nov 

27.2 

60.0 

18.8 

8.1 

Dec 

20 .8 

46.0 

13.2 

7.9 

Table  A9.  Eudora,  second  cel  1,  44  days. 

FI  nal  effluent 

No 

Input,  total 

Temp 

» 

total  N  <mq/t) 

N  <  mq/t) 

CC) 

pH 

Jan 

16.7 

25.1 

4.5 

8.1 

Feb 

16.5 

21,8 

3.9 

8.3 

Nter 

15.6 

24.3 

4.5 

8.5 

Apr 

18.8 

24.3 

7.0 

8.6 

’ 

May 

18.0 

24.8 

16.3 

8.6 

Jun 

9.8 

30.6 

22.0 

8.5 

t 

Jul 

15.9 

28.0 

24.6 

8.8 

Aug 

18.3 

16.4 

26.7 

8.8 

Sep 

9.1 

24.1 

24.6 

8.5 

Oct 

14.0 

21.6 

20.3 

8.5 

Nov 

15.5 

17.8 

14.6 

8.2 

■»n=i 

13.3 

23.4 

8.1 

8.0 

1 

21 

1 

Table  A10.  Eudora,  first  and  second  cells,  137  days, 


Final  effluent 
total  N  (mqA) 

No 

Input,  total 

N  ( mg  A) 

Tamp 

CC) 

pH 

Jan 

16.7 

53.0 

12.5 

8.1 

Feb 

16.5 

45.0 

8.2 

8.0 

Msr 

15.6 

47.5 

5.5 

8.0 

Apr 

18.8 

49.0 

4.1 

8.2 

May 

18.0 

48.5 

8.6 

8.2 

Jun 

9.8 

52.0 

12.7 

8.3 

J  ul 

15.9 

52.0 

17.6 

8.4 

Aug 

18.3 

59.0 

21.8 

8.4 

Sep 

9.1 

55.5 

24.0 

8.4 

Oct 

14.0 

39.0 

23.7 

8.5 

Nov 

15.5 

48.0 

21.3 

8.3 

Dec 

13.3 

60.5 

17.0 

8.2 

Table  All.  Corlnne,  three  cel  Is,  42  days. 


Final  effluent 
total  N  (mqA) 

Ho 

Input,  total 

N  ( mq  A ) 

Temp 

CC) 

pH 

Alka 1 1 nlty 
(CaC03) 

(  mg  A) 

Jan 

7.2 

23.4 

1.9 

9.3 

530 

Feb 

10.6 

19.6 

1.2 

9.0 

550 

Mir 

11.9 

14.1 

2.8 

B.8 

530 

Apr 

8.7 

15.3 

7.7 

9.0 

500 

►toy 

5.4 

16.1 

10.0 

9.5 

510 

Jun 

5.0 

16.1 

14.6 

9.6 

550 

J  ul 

4.1 

11.2 

20.0 

9.4 

590 

Aug 

5.9 

8.2 

21.1 

9.4 

600 

Sep 

7.0 

10.3 

17.8 

9.5 

620 

Oct 

7.2 

14.6 

12.2 

9.5 

610 

Nov 

6.8 

14.6 

5.7 

9.6 

550 

Dec 

4.6 

17.6 

2.6 

9.6 

540 

Table  A12.  Corlnne,  first  cell,  27  days. 


Final  effluent 
total  N  (mg A) 

No 

Input,  total 

N  (mqA) 

Temp 

CC) 

pH 

Jan 

14.0 

23.4 

2.1 

8.8 

Feb 

15.4 

19.6 

1.6 

8.4 

Ifcr 

12.1 

14.1 

3.3 

8.1 

Apr 

9.4 

15.3 

6.9 

8.5 

May 

7.5 

16.1 

10.1 

9.3 

Jun 

6.9 

16.1 

14.7 

9.5 

J  ul 

6.7 

11.2 

15.0 

9.4 

Aug 

8.7 

8.2 

20.7 

9.4 

Sep 

8.2 

10.3 

17.6 

9.4 

Oct 

8.4 

14.6 

12.2 

9.4 

Nov 

8.5 

14.6 

5.7 

9.5 

Dec 

8.0 

17.6 

2.8 

9.5 

22 


Table  A13.  Corlnne,  first  and  second  cells,  34.5  days, 


Final  effluent 
total  N  (mqA.) 

No 

Input,  total 

N  (mq/O 

Temp 

CC> 

pH 

Jan 

8.9 

23.4 

1.9 

9.2 

Feb 

12.3 

19.6 

1.3 

8.8 

ftor 

12.4 

14.1 

3.0 

8.8 

Apr 

8.4 

15.3 

6.2 

9.1 

May 

6.1 

16.1 

10.0 

9.4 

Jtin 

4.9 

16.1 

14.8 

9.6 

Jul 

5.0 

11.2 

20.2 

9.4 

Aug 

7.6 

8.2 

21.0 

9.4 

Sep 

8.2 

10.3 

17.8 

9.5 

Oct 

7.6 

14.6 

12.3 

9.5 

Nov 

6.9 

14.6 

5.8 

9.6 

Dec 

5.8 

17.6 

2.7 

9.6 

23 
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APPENDIX  B:  DATA  USED  FOR  MODEL  VALIDATION 


5 

Table  B1.  Gulfport,  one  cell,  5  days. 


Date 

Final  effluent 
total  N  (mg/L) 

No 

Input,  total 

N  ( mq  A.) 

Temp 

CC) 

pH 

Dec  75 

14.4 

15.4 

10.9 

7.2 

Jan  76 

15.8 

15.5 

12.0 

7.2 

Feb 

15.0 

13.5 

18.8 

7.1 

Mar 

12.5 

16.3 

15.3 

7.0 

Apr 

15.0 

16.7 

21.4 

7.2 

Hay 

14.0 

17.6 

22.4 

7.6 

Jun 

10.0 

11.5 

26.9 

7.3 

Jul 

9.7 

10.8 

28.0 

7.3 

Aug 

9.9 

10.0 

28.1 

7.2 

Sep 

9.4 

9.7 

25.0 

7.2 

Oct 

13.8 

13.1 

23.6 

6.9 

Nov 

11.6 

12.4 

18.0 

7.1 

Table  ^2.  Drayton  Valley,  final  two  cells,  69  days  (5.5  ha,  1.52  m 
deep). 


Date 

Final  effluent 
NHu+  (mq/L) 

No 

Input,  total 
NHu  (mq/L) 

T  emp 
(*C> 

pH 

A Ikal Ini ty 
(mq/L) 

Ice 

depth 

(m) 

Jul  62 

5.0 

20.6 

17.0 

8.5 

560 

0 

Aug 

3.5 

14.3 

20.7 

9.2 

435 

0 

Nov 

6.3 

13.5 

11.3 

9.0 

640 

0.05 

Jan  63 

35.0 

47.3 

1.2 

8.3 

930 

0.76 

Mar 

33.0 

47.8 

1.3 

8.0 

415 

0.61 

Apr 

26.0 

31.0 

1.8 

7.7 

500 

0.61 

May 

12.0 

28.0 

6.9 

7.9 

465 

0 

Dec 

18.0 

32.0 

1.3 

8.2 

790 

0.15 

Mar  64 

42.0 

48.0 

1.5 

7.9 

855 

0.60 

Table  B3^  Stettler,  first  long  detention  cell,  75  days  (5.34  ha,  1.52 
m  deep). 


Date* 

Final  effluent 
ammonia 
(mq/L) 

1  nput 
ammonia 
(mq/l) 

Temp 

CC) 

P_H 

Alka 1 1 nlty 
( mq/L) 

Ice 

depth 

(m) 

J  ul 

17 

49. 5t 

20.0 

8.7 

612 

0 

Aug 

14 

49. 5t 

19.0 

8.9 

835 

0 

Nov 

37 

49. 5t 

5.0 

8.5 

870 

0 

Dec 

49 

49. 5t 

1.5 

7.9 

1115 

0.30 

Jan 

47 

49. 5t 

0.7 

7.7 

1277 

0.91 

Feb 

55 

49. 5t 

1.0 

8.0 

1345 

0.61 

Mir 

58 

49. 6t 

1.0 

7.6 

1030 

1.1 

May 

28 

49. 5t 

12.0 

7.7 

950 

0 

Jun 

18 

49. 5t 

18.5 

8.1 

965 

0 

•Values  are  averages  tor  month  Indicated  tor  1961-64.  Data  not  col- 
lected  In  other  months, 
tAverage  of  18  samples,  1961-64. 
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Table  B4. 

Ini  stall,  first 

long  detention  cell 

,  22. 

6  cays  (1.18 

ha,  1. 

m  deep). 

Final  effluent 

Input 

Ice 

ammonia 

ammonia 

Temp 

Alka 1 Inity 

depth 

Date* 

KE  EJM 

CO 

pH 

■a  h— 

( m) 

Dec 

39 

47 

0.8 

7.6 

488 

0.46 

Jan 

54 

46 

1.0 

7.5 

512 

0.66 

Feb 

46 

40 

1.1 

7.4 

525 

0.76 

Mar 

49 

40 

1.8 

7.4 

508 

0.91 

Apr 

36 

40 

2.9 

7.4 

330 

0 

May 

34 

39 

8.3 

7.5 

420 

0 

Jul 

34 

37 

18.0 

8.2 

348 

0 

Sep 

24 

34 

18.0 

8.1 

380 

0 

Nov 

25 

43 

12.0 

7.9 

415 

0.03 

Average  42.2 

•Values  are  averages  for  month  Indicated  for  1961-64.  Data  not  col¬ 
lected  In  other  months. 


Table  B5.  Laeombe,  one  cel  I  (6.73  ha.  1.52  m  deep).  Drained  In  late 
October,  refill  and  storage  October-torch;  discharge  equal  to  Input 
torch  to  September. 


a 

« 


Date* 

Final  effluent 
ammonia 
( mq  A) 

1  nput 
ammonia 
(mqA) 

Tenp 

CC) 

pH 

Alka 1 1 nlty 
( mq  A) 

Ice 

depth 

(m) 

Jan 

83 

51.25t 

0.73 

8.2 

1775 

0.46 

Feb 

57 

51.25t 

1.0 

7.8 

882 

0.88 

tor 

39 

51.25t 

1.0 

7.9 

963 

0.91 

Apr 

28 

51.25t 

5.0 

7.4 

555 

0 

May 

12 

5l.25t 

24.0 

9.4 

600 

0 

Jul 

15 

51.25t 

18.0 

8.3 

535 

0 

Aug 

10 

51 ,25t 

17.0 

9.3 

490 

0 

Nov 

22 

51.25t 

4.0 

8.1 

760 

0.06 

Dec 

47 

51 ,25t 

1.0 

7.9 

1200 

0. 30 

•Values  are  averages  for  month  Indicated  for  l'56l-64.  Data  not  col- 
lected  In  other  months. 

tAverage  of  four  samples  of  raw  sewage  taken  In  1963-64. 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 


Reed,  Sherwood  C. 

Nitrogen  removal  in  wastewater  ponds  /  by  Sherwood 
C.  Reed.  Hanover,  N.H.:  Cold  Regions  Research  and 
Engineering  Laboratory;  Springfield,  Va. :  available 
from  National  Technical  Information  Service,  1984. 

v,  33  p.,  illus.;  28  cm.  (  CRREL  Report  84-13.  ) 

Prepared  for  Office  of  the  Chief  of  Engineers 
by  Corps  of  Engineers,  U.S.  Army  Cold  Regions  Re¬ 
search  and  Engineering  Laboratory  under  Da  Project 
4A762730AT42. 

Bibliography:  p.  17. 
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5.  Waste  water.  I.  United  States.  Army.  Corps  of 
Engineers.  II.  Cold  Regions  Research  and  Engineer¬ 
ing  Laboratory,  Hanover,  N.H.  III.  Series:  CRREL 
Report  84-13. 


